Genome-Wide Analysis of Differentially Expressed miRNAs and Their Associated Regulatory Networks in Lenses Deficient for the Congenital Cataract-Linked Tudor Domain Containing Protein TDRD7.

Mutations/deficiency of TDRD7, encoding a tudor domain protein involved in post-transcriptional gene expression control, causes early onset cataract in humans. While Tdrd7 is implicated in the control of key lens mRNAs, the impact of Tdrd7 deficiency on microRNAs (miRNAs) and how this contributes to transcriptome misexpression and to cataracts, is undefined. We address this critical knowledge-gap by investigating Tdrd7-targeted knockout (Tdrd7-/-) mice that exhibit fully penetrant juvenile cataracts. We performed Affymetrix miRNA 3.0 microarray analysis on Tdrd7-/- mouse lenses at postnatal day (P) 4, a stage preceding cataract formation. This analysis identifies 22 miRNAs [14 over-expressed (miR-15a, miR-19a, miR-138, miR-328, miR-339, miR-345, miR-378b, miR-384, miR-467a, miR-1224, miR-1935, miR-1946a, miR-3102, miR-3107), 8 reduced (let-7b, miR-34c, miR-298, miR-382, miR-409, miR-1198, miR-1947, miR-3092)] to be significantly misexpressed (fold-change ≥ ± 1.2, p-value < 0.05) in Tdrd7-/- lenses. To understand how these misexpressed miRNAs impact Tdrd7-/- cataract, we predicted their mRNA targets and examined their misexpression upon Tdrd7-deficiency by performing comparative transcriptomics analysis on P4 and P30 Tdrd7-/- lens. To prioritize these target mRNAs, we used various stringency filters (e.g., fold-change in Tdrd7-/- lens, iSyTE-based lens-enriched expression) and identified 98 reduced and 89 elevated mRNA targets for overexpressed and reduced miRNAs, respectively, which were classified as "top-priority" "high-priority," and "promising" candidates. For Tdrd7-/- lens overexpressed miRNAs, this approach identified 18 top-priority reduced target mRNAs: Alad, Ankrd46, Ceacam10, Dgat2, Ednrb, H2-Eb1, Klhl22, Lin7a, Loxl1, Lpin1, Npc1, Olfm1, Ppm1e, Ppp1r1a, Rgs8, Shisa4, Snx22 and Wnk2. Majority of these targets were also altered in other gene-specific perturbation mouse models (e.g., Brg1, E2f1/E2f2/E2f3, Foxe3, Hsf4, Klf4, Mafg/Mafk, Notch) of lens defects/cataract, suggesting their importance to lens biology. Gene ontology (GO) provided further insight into their relevance to lens pathology. For example, the Tdrd7-deficient lens capsule defect may be explained by reduced mRNA targets (e.g., Col4a3, Loxl1, Timp2, Timp3) associated with "basement membrane". GO analysis also identified new genes (e.g., Casz1, Rasgrp1) recently linked to lens biology/pathology. Together, these analyses define a new Tdrd7-downstream miRNA-mRNA network, in turn, uncovering several new mRNA targets and their associated pathways relevant to lens biology and offering molecular insights into the pathology of congenital cataract.

The Tudor-domain protein TDRD7, mutated in congenital cataract, controls the heat shock protein HSPB1 (HSP27) and lens fiber cell morphology.

Mutations of the RNA granule component TDRD7 (OMIM: 611258) cause pediatric cataract. We applied an integrated approach to uncover the molecular pathology of cataract in Tdrd7-/- mice. Early postnatal Tdrd7-/- animals precipitously develop cataract suggesting a global-level breakdown/misregulation of key cellular processes. High-throughput RNA sequencing integrated with iSyTE-bioinformatics analysis identified the molecular chaperone and cytoskeletal modulator, HSPB1, among high-priority downregulated candidates in Tdrd7-/- lens. A protein fluorescence two-dimensional difference in-gel electrophoresis (2D-DIGE)-coupled mass spectrometry screen also identified HSPB1 downregulation, offering independent support for its importance to Tdrd7-/- cataractogenesis. Lens fiber cells normally undergo nuclear degradation for transparency, posing a challenge: how is their cell morphology, also critical for transparency, controlled post-nuclear degradation? HSPB1 functions in cytoskeletal maintenance, and its reduction in Tdrd7-/- lens precedes cataract, suggesting cytoskeletal defects may contribute to Tdrd7-/- cataract. In agreement, scanning electron microscopy (SEM) revealed abnormal fiber cell morphology in Tdrd7-/- lenses. Further, abnormal phalloidin and wheat germ agglutinin (WGA) staining of Tdrd7-/- fiber cells, particularly those exhibiting nuclear degradation, reveals distinct regulatory mechanisms control F-actin cytoskeletal and/or membrane maintenance in post-organelle degradation maturation stage fiber cells. Indeed, RNA immunoprecipitation identified Hspb1 mRNA in wild-type lens lysate TDRD7-pulldowns, and single-molecule RNA imaging showed co-localization of TDRD7 protein with cytoplasmic Hspb1 mRNA in differentiating fiber cells, suggesting that TDRD7-ribonucleoprotein complexes may be involved in optimal buildup of key factors. Finally, Hspb1 knockdown in Xenopus causes eye/lens defects. Together, these data uncover TDRD7's novel upstream role in elevation of stress-responsive chaperones for cytoskeletal maintenance in post-nuclear degradation lens fiber cells, perturbation of which causes early-onset cataracts.

Population genetic correlates of declining transmission in a human pathogen.

Pathogen control programs provide a valuable, but rarely exploited, opportunity to directly examine the relationship between population decline and population genetics. We investigated the impact of an ~12-fold decline in transmission on the population genetics of Plasmodium falciparum infections (n = 1731) sampled from four clinics on the Thai-Burma border over 10 years and genotyped using 96 genome-wide SNPs. The most striking associated genetic change was a reduction in the frequency of infections containing multiple parasite genotypes from 63% in 2001 to 14% in 2010 (P = 3 × 10(-15)). Two measures of the clonal composition of populations (genotypic richness and the ß-parameter of the Pareto distribution) declined over time as more people were infected by parasites with identical multilocus genotypes, consistent with increased selfing and a reduction in the rate at which multilocus genotypes are broken apart by recombination. We predicted that the reduction in transmission, multiple clone carriage and outbreeding would be mirrored by an increased influence of genetic drift. However, geographical differentiation and expected heterozygosity remained stable across the sampling period. Furthermore, N(e) estimates derived from allele frequencies fluctuation between years remained high (582 to ∞) and showed no downward trend. These results demonstrate how genetic data can compliment epidemiological assessments of infectious disease control programs. The temporal changes in a single declining population parallel to those seen in comparisons of parasite genetics in regions of differing endemicity, strongly supporting the notion that reduced opportunity for outbreeding is the key driver of these patterns.

A major genome region underlying artemisinin resistance in malaria.

Evolving resistance to artemisinin-based compounds threatens to derail attempts to control malaria. Resistance has been confirmed in western Cambodia and has recently emerged in western Thailand, but is absent from neighboring Laos. Artemisinin resistance results in reduced parasite clearance rates (CRs) after treatment. We used a two-phase strategy to identify genome region(s) underlying this ongoing selective event. Geographical differentiation and haplotype structure at 6969 polymorphic single-nucleotide polymorphisms (SNPs) in 91 parasites from Cambodia, Thailand, and Laos identified 33 genome regions under strong selection. We screened SNPs and microsatellites within these regions in 715 parasites from Thailand, identifying a selective sweep on chromosome 13 that shows strong association (P = 10(-6) to 10(-12)) with slow CRs, illustrating the efficacy of targeted association for identifying the genetic basis of adaptive traits.

Emergence of artemisinin-resistant malaria on the western border of Thailand: a longitudinal study.

BACKGROUND: Artemisinin-resistant falciparum malaria has arisen in western Cambodia. A concerted international effort is underway to contain artemisinin-resistant Plasmodium falciparum, but containment strategies are dependent on whether resistance has emerged elsewhere. We aimed to establish whether artemisinin resistance has spread or emerged on the Thailand-Myanmar (Burma) border. METHODS: In malaria clinics located along the northwestern border of Thailand, we measured six hourly parasite counts in patients with uncomplicated hyperparasitaemic falciparum malaria (≥4% infected red blood cells) who had been given various oral artesunate-containing regimens since 2001. Parasite clearance half-lives were estimated and parasites were genotyped for 93 single nucleotide polymorphisms. FINDINGS: 3202 patients were studied between 2001 and 2010. Parasite clearance half-lives lengthened from a geometric mean of 2·6 h (95% CI 2·5-2·7) in 2001, to 3·7 h (3·6-3·8) in 2010, compared with a mean of 5·5 h (5·2-5·9) in 119 patients in western Cambodia measured between 2007 and 2010. The proportion of slow-clearing infections (half-life ≥6·2 h) increased from 0·6% in 2001, to 20% in 2010, compared with 42% in western Cambodia between 2007 and 2010. Of 1583 infections genotyped, 148 multilocus parasite genotypes were identified, each of which infected between two and 13 patients. The proportion of variation in parasite clearance attributable to parasite genetics increased from 30% between 2001 and 2004, to 66% between 2007 and 2010. INTERPRETATION: Genetically determined artemisinin resistance in P falciparum emerged along the Thailand-Myanmar border at least 8 years ago and has since increased substantially. At this rate of increase, resistance will reach rates reported in western Cambodia in 2-6 years. FUNDING: The Wellcome Trust and National Institutes of Health.

Transmission and cross-mating of high-level resistance Plasmodium falciparum dihydrofolate reductase haplotypes in The Gambia.

A high-level pyrimethamine resistance Plasmodium falciparum lineage with triple dihydrofolate reductase (dhfr) mutations prevails across Africa. However, additional minority lineages were seen. We examined transmission success of mutant dhfr haplotypes among 22 children in The Gambia and 60 infected Anopheles gambiae mosquitoes fed on their blood. Additional polymorphic genes of the gametocyte-specific protein (pfg377) and merozoite surface protein-1 (MSP-1) were examined. Similarities were seen between pfg377 and MSP-1 alleles in children and mosquitoes and evidence of cross-mating between different parasite genotypes was seen in some infected mosquitoes, reflecting high transmission success of existing clones. With regard to dhfr, 16 haplotypes were seen among the children: 2 carried double mutations and 14 carried triple mutations. However, only nine haplotypes, all with triple mutations, were detected among mosquitoes. A single triple-mutant dhfr haplotype, similar to that in other countries in Africa, predominated among children (42%) and mosquitoes (60%), supporting the hypothesis of migration of this haplotype across Africa. However, evidence of cross-mating between the above haplotypes signifies the role of local evolution.

Distinct haplotypes of dhfr and dhps among Plasmodium falciparum isolates in an area of high level of sulfadoxine-pyrimethamine (SP) resistance in eastern Sudan.

Typing of polymorphic microsatellites that are linked to drug resistance genes has shed light on the origin and pattern of spread of some anti-malarial drugs. Recent surveys revealed spread of a high-level pyrimethemine resistant lineage of Plasmodium falciparum, of Asian origin, across Africa. Here, we examined mutations in dihydrofolate reductase, dhfr [chromsosome 4], the dihydropteroate synthase, dhps [chromosome 8] associated with resistance to sulfadoxine-pyrimethamine (SP), and neighboring microsatellites among P. falciparum isolates in Asar village, eastern Sudan. This area lies at the fringes of malaria endemicity, where the remote P. falciparum parasites have some distinct genetic characteristics. Overall, 89% (84/94) of the examined isolates carried double mutations at dhfr (N51I and S108N), but the 59R and I164L mutations were not seen. Similarly, the majority, 43% (35/81) of the isolates carried double mutations at dhps (437G, 540E). Analysis of neighboring microsatellites revealed one major dhfr haplotype with mutations (51I, 108N) and one dhps haplotype with mutations (436S, 437G, 540E). These haplotypes differ from the major ones thought to drive resistance to SP across Africa. The resistant haplotypes of dhfr and dhps, in Asar, share some microsatellites with the wild genotypes suggesting that they were generated locally. Among isolates successfully examined, 40% shared identical haplotypes of the 2 loci, comprising a dominant resistant lineage. Undoubtedly, this lineage plays an important role in clinical failure to SP in this area.